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ABSTRACT 

Gravitational flexion has been introduced as a technique by which one can map out 
and study substructure in clusters of galaxies. Previous analyses involving flexion 
have measured the individual galaxy-galaxy flexion signal, or used either parametric 
techniques or a KSB-type inversion to reconstruct the mass distribution in Abell 1689. 
In this paper, we present an aperture mass statistic for flexion, and apply it to the 
lensed images of background galaxies obtained by ray-tracing simulations through a 
simple analytic mass distribution and through a galaxy cluster from the Millennium 
simulation. We show that this method is effective at detecting and accurately tracing 
structure within clusters of galaxies on sub-arcminute scales with high signal-to-noise 
even using a moderate background source number density and image resolution. In 
addition, the method provides much more information about both the overall shape 
and the small-scale structure of a cluster of galaxies than can be achieved through a 
weak lensing mass reconstruction using gravitational shear data. Lastly, we discuss how 
the zero-points of the aperture mass might be used to infer the masses of structures 
identified using this method. 

Key words: cosmology:observations - cosmology: dark matter - galax- 
ies: clusters: general - gravitational lensing 



1 INTRODUCTION 

Galaxy clusters are the most massive bound structures in the 
universe, some in excess of IQ^^Mq, with about 90% of their 
mass in the form of dark matter. The mass function of clus- 
ters provides a sensitive test of our cosmological model (e.g. 
Bahcall & Fan 1998; Eke et al. 1998), and the mass func- 
tion of substructure provides important insight into galaxy 
formation in dense environments (e.g. De Lucia et al. 2004; 
Taylor & Babul 2005). 

In a universe dominated by Cold Dark Matter (CDM) 
such as our concordance cosmology, ACDM, hierarchical 
structure formation takes place with the growth of collapsed 
objects progressing via merging of smaller objects. In non- 
hierarchical structure formation, taking place in a universe 
dominated by Hot Dark Matter (HDM) for example, the first 
haloes form via monolithic collapse. Recent N-body simula- 
tions of halo formation in a HDM dominated universe by 
Wang & White (2008) show that it is the mass substructure 
content of haloes - rather than most other halo properties 
such as shape or spin parameter - that is markedly different 
from that seen in CDM. 

Gravitational lensing, being insensitive to whether mat- 
ter is luminous or dark and to its dynamical state, is an 
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ideal probe of substructure. Natarajan, De Lucia & Springel 

(2007) have obtained constraints on the substructure in five 
massive galaxy clusters using weak and strong gravitational 
lensing observations, finding comparable levels (~ 10 — 20%) 
to that seen in high-resolution simulated clusters. 

In addition, Leonard et al. (2007) and Okura et al. 

(2008) have used flexion measurements to detect substruc- 
ture in Abell 1689. While the Okura et al. reconstruction 
uses a nonparametric, KSB-type (Kaiser, Squires & Broad- 
hurst, 1995) reconstruction technique, the Leonard et al. re- 
construction relies on a parametric modelling of known clus- 
ter members. Each of these methods has its disadvantages; 
using the finite field inversion technique of Okura et al., it is 
difficult to accurately characterise the noise properties of the 
mass maps produced, and thus to assess the significance of 
the detections, whilst the Leonard et al. reconstruction re- 
quires detailed knowledge of the locations of the structures 
responsible for the lensing. Since flexion studies in galaxy 
clusters are highly sensitive to the substructure within the 
cluster, these studies cannot be compared to weak or strong 
lensing reconstructions, or to mass models of the cluster as 
a whole, to assess their performance. 

In this paper, we derive an aperture mass statistic for 
gravitational flexion in direct analogy to that for gravita- 
tional shear (e.g. Schneider 1996; Schneider et al. 1998), to 
enable us to use measurements of flexion to determine the lo- 
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cations and statistical significance of mass peaks in clusters 
of galaxies without the need for any parametric modelling or 
comparison with complementary lensing studies. We begin 
with a review of the flexion formalism, outline how flexion 
is measured and techniques that have previously been used 
to map the mass distribution in Abell 1689. In Section 3 the 
aperture statistic for flexion is derived, and in Section 4 we 
describe the ray-tracing simulations that have been carried 
out and how the aperture statistic is calculated from the syn- 
thetic data. Section 5 presents results on the performance 
of the statistic, and in Section 6 we discuss our findings and 
conclude. 

Throughout this paper wc use a matter density param- 
eter Om = 0.27, dark energy density parameter Qa ~ 0.73 
(with equation of state parameter w = —1) and Hubble pa- 
rameter Ho = 71km s~^ Mpc~^. 



2 REVIEW OF FLEXION FORMALISM 

In traditional weak lensing studies, the lens equation is ap- 
proximated as linear, and lensed images exhibit a purely 

elliptical distortion aligned tangential to the lens. Linearisa- 
tion of the lens equation assumes that there is no variation 
of the lens field over the scale of the lensed image. If this 
field is allowed to vary smoothly, the lens equation becomes 
non-linear (see Goldberg & Bacon 2005): 

Pi ~ AijOj + -DijkdjOk, (1) 

where (3 is the coordinate in the source plane, is the lensed 
coordinate, A is the magnification matrix 

K is the convergence, 71 , 72 are the components of the (com- 
plex) gravitational shear (a polar quantity), and Dijk = 
dkAij. The D operators can be related to two measurable 
quantities, first and second fiexion, by 

2 I, ^2+02 J^i-Gi J' 

^^^^--2 I, ^1-01 3^2-02 )' 

where = J^i + ij^2 = On, Q = Qi + iQ2 = dj, and d = 
di + 182 is the differential operator defined in Bacon et al. 

(2006) . 

First flexion, .F, transforms as a vector, directly probing 
the gradient of the convergence, and gives rise to a skcwness 
in the light distribution of the lensed image aligned radially 
with respect to the lens. Second fiexion, Q has m = 3 rota- 
tional symmetry and gives rise to an arciness in the lensed 
image, aligned tangential to the lens. Goldberg & Leonard 

(2007) and Leonard et al. (2007) found that second flexion 
is signiflcantly more difficult to measure than first flexion, 
thus for the purposes of this paper, we will consider only 
first flexion. 

There are two distinct methods by which flexion can be 
measured. The flrst, described in detail in Goldberg & Bacon 
(2005), Bacon et al. (2006) and Massey et al. (2007) involves 
the decomposition of lensed images into an orthogonal basis 



set, shapelets (see e.g. Bernstein & Jarvis 2002, Refregier 
2003, Refregier & Bacon 2003, Massey & Refregier 2005), 
which are related to two-dimensional reduced Hermite or 
Laguerre polynomials. One advantage of this method is that 
the various lensing operators can be expressed rather simply 
in terms of the quantum-mechanical raising and lowering op- 
erators, and produce rather compact transfers of power be- 
tween shapelet modes. Another advantage is that an explicit 
deconvolution of the telescope point spread function (PSF) 
can be carried out rather straightforwardly in this method. 

An alternative technique has been formulated by Okura 
et al. (2007) and subsequently refined and extended in 
Goldberg & Leonard (2007) and Okura et al. (2008). This 
technique involves measuring higher order moments of the 
brightness distribution of the lensed image (HOLIGs), and 
using these moments to estimate the flexion. While less com- 
pact on paper, the HOLIGs method has the advantage that 
the computational anatysis time is dramatically shorter than 
that involved with shapelets, particularly on large, well- 
resolved lensed images. In addition, incorporating a weight- 
ing function in the measurement of the HOLIGs reduces the 
impact of noise in the image, as well as light contamination 
from nearby sources, and the latest work by Okura et al. 
(2008) describes a method whereby the effects of both an 
isotropic and an anisotropic PSF can be corrected for. 

The shapelets technique has been used successfully to 
measure galaxy-galaxy flexion in a sample of field galax- 
ies in the Deep Lens Survey (Goldberg and Bacon 2005), 
while two different implementations of the HOLIGs formal- 
ism have been used to measure flexion in images of Abell 
1689 (Leonard et al. 2007, Okura et al. 2008), where it has 
been shown that flexion is particularly sensitive to substruc- 
ture in clusters of galaxies. While the Leonard et al. recon- 
struction includes a description of the noise properties of 
their reconstructed convergence map, the technique relies 
on accurate knowledge of the locations of the cluster mem- 
bers responsible for the measured flexion signal, and takes 
no account of the smooth component of the cluster poten- 
tial, nor does it account for the possible presence of dark 
haloes within the cluster. The Okura et al. reconstruction 
has the advantage of using a nonparamctric technique, and 
thus not requiring any a priori assumptions about the clus- 
ter. However, it is very difficult with this technique to accu- 
rately describe the noise properties of the convergence map 
generated. 

It is in this context that an aperture mass statistic, 
which includes a straightforward description of the noise 
properties of the mass peaks detected, is important. 



3 THE APERTURE MASS STATISTIC 

To derive the aperture mass statistic for flexion, we closely 
follow the work of Schneider (1996), in which the formalism 
for the generalised aperture mass statistic for shear was first 
laid out. The aperture mass is defined as: 

m{xo) — J d^x k(x + xo) '(u(|x|). (4) 

We aim to express this statistic in terms of some measurable 
quantity, namely the measured fiexion. 

The convergence is related to the first fiexion by (see 
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xi + ix2- Thus, we can redefine the aperture mass 



j (fx wix) J d^x' £;>(x-x'+xo):'^(x'; 



/ 



+ / d X ■w{x)ko 



(7) 



Schneider (1996) notes that the aperture mass can be made 
independent of the constant quantity kq if we require that 
the mass filter function w{x) is compensated, i.e. 



X w{x)dx — 0. 



(8) 



Under this condition, the aperture mass for fiexion can 
be redefined in terms of the "E-mode" (radially afigned) 
fiexion as follows (see Appendix |B] for the complete deriva- 
tion): 
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In direct analogy with the shear aperture mass statis- 
tic, the expected signal to noise ratio, S, achievable for a 
given fiexion filter function can be calculated by taking an 
ensemble average over the probability distribution for the 
background galaxy positions, which gives 



-{x) {Te) {x) dx 
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where n is the number density of background sources within 
the aperture, and ctjf is the dispersion in fiexion measure- 
ments. Goldberg & Leonard (2007) found <ya\T\ = 0.03, im- 
plying ajr ~ 0.28/". We note that this is likely to be an over- 
estimate of the true, unlensed dispersion in fiexion values, 
however, as this estimate is based on the fiexion measured 
in galaxies within a moderate lensing field. 

3.1 Choosing a Filter 

We consider two difi'erent families of filter functions for the 
fiexion aperture mass statistic. These filters are able to at- 
tain a very similar peak signal-to-noise; however, the depen- 
dence of the signal to noise ratio on various factors (such as 
the scale of the lens substructures and the size of the aper- 
ture being used) differs significantly between the two sets of 
filters, thus their respective domains of applicability differ. 

It is important to note that the signal to noise attain- 
able for a given lens profile is maximised when the flexion 
filter function traces the expected flexion signal. Thus the 



optimal choice of fllter function is strongly dependent on 
the flexion profiles of the structures being studied. The fil- 
ter functions presented in this paper have not been chosen 
to be optimal for a given lens profile; rather they are de- 
signed to have rather broad applicability, thus allowing us 
to detect structures without requiring a priori knowledge of 
the profiles of these structures. 



3.1.1 Piecewise- continuous Filters 

Schneider (1996) describes a generalised piecewise- 
continuous, compensated mass filter function, from which 
he derives a shear filter function: 

1 xKviR 
W.,,(x) = < I , '^'-"^ — C I :^iR<x<i^2R , 



[^{R-x)\x-aR) 



(12) 

where v\R is an inner radius, usually taken to be small com- 
pared to the size of the aperture, and V2R is an outer radius, 
usually taken to be close to the aperture scale R. The con- 
stants a, b and c are calculated using the constraints that 
w^{x) and dw.y(x) / dx must be continuous at x — uiR and 
X = V2R, and that w^(x) must be compensated. 

Noting that the mass filter function for flexion is roughly 
proportional to the derivative of the mass fllter function for 
shear, we choose 



1 



((a:-,.lfl)2 + (i,iii)2)* 
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C viR<x<U2R , 



V2R<x<R, 



(13) 

where, again, the constants can be computed by requiring 
that w{x) and dw{x)/dx be continuous at a: = uiR and 
X = i'2R, and that w{x) be compensated. Applying equation 
fTol we flnd that 



2' 



Qr^viR < a; ^ V2R) 
{u^R)^ f 3 + 
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-20W-X (4(^' - (^2^)') - 5(2 + a)R{x^ - {v2Rf)) 
(2(1 + 2a)(x3 - {u2Rf) - 3aR{x^ - {v2Rf)) . 

(14) 

Fig.[l]shows the mass and flexion filter functions derived 
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Figure 1. Left Panel: The piecewise-continuous mass filter func- 
tion described in § |3.1.1| plotted as a function of x/ij, wliere R is 
the aperture radius, for 1^2 = 0.95. Right Panel: The correspond- 
ing flexion filter functions. 



Figure 2. Left Panel: The polynomial mass filter function de- 
scribed in § |3.1.2| plotted as a function oi x/R. Right Panel: The 
corresponding flexion filter functions. 



above for various choices of ui, taking 1/2 = 0.95. The flexion 
filter drops of roughly as for uiR ^ x ^ z^2-R, thus is 
not optimised for either an isothermal profile or an NFW 
profile (and, in fact, would only be optimal for a profile with 
K,{x) oc log(a;)). However, as it has a rather wide distribution, 
we expect to be able to detect structures with a broad range 
of profiles using this filter. 



3.1.2 Polynomial Filters 

A more compact family of (continuous) mass filter functions 
was described by Schneider et al. (1998) as: 
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(15) 



where ^ is a factor arising from the normalisation of the 
associated fiexion filter function. This family of functions 
drops off with order I as x ^ R, and were designed to be 
applicable for cosmic shear studies. Thus, the filter functions 
are not optimised, but should be sensitive to a range of mass 
profiles. 

The associated flexion filter functions are given by 
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It is clear from dimensionality considerations that Qj^ must 
have dimensions of length, as the product J-Qr must be 
dimensionless and flexion itself carries dimensions of inverse 
length. Thus, we choose the normalisation constant. A, such 
that 



2-K 



Jo 
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Thus, 
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These mass and flexion filter functions are plotted in 
Fig. [2] for various values of I. 



3.2 Expected Signal to Noise 

The expected signal to noise achievable for a given filter 
function can be calculated according to equation |11[ This 
turns out to be quite a straightforward calculation using the 
polynomial flexion filter functions (equation |16[ ) and under 
the assumption of a singular isothermal sphere lens proflle 



{J-e{x) = 9e/2x'^), where we find the signal to noise to be 
given by: 

where n is the background source density, ajr is the intrinsic 
fiexion dispersion, and Oe is the Einstein radius of the lens, 
given by 



(20) 



where Ds is the angular diameter distance to the source 
plane and Dis is the angular diameter distance between the 
lens and source planes. 

Thus, we can see that the expected signal to noise for 
this set of filters decreases with increasing aperture size. This 
would seem to imply that, in order to optimise the signal to 
noise, the apertures should be made as small as possible. 
However, decreasing the aperture size decreases the number 
of sources within each aperture. In a noisy image where the 
fiexion is not measured perfectly, increasing the number of 
sources within an aperture will increase the statistical signif- 
icance of the measurement. In other words, having a small 
number of (noisy) flexion measurements within an aperture 
will artiflcially increase the value of ajr within that aper- 
ture, thus decreasing the signal to noise achievable within 
the aperture. 

In addition, a small aperture size will inevitably result 
in many apertures over the field having no sources in them, 
particularly when using data with a low background source 
count. This will result in a patchy measured aperture mass 
with discontinuities seen where there are apertures contain- 
ing no sources. It is thus necessary to use larger apertures 
than the signal to noise calculation might imply are optimal. 

The above calculation also assumes that the fiexion sig- 
nal is measurable at a; = 0, and indeed more heavily weights 
the signal close to the centre of the aperture. Now, in prac- 
tice, we will be unable to measure the flexion at very small 
distances from the centre of the lens because the lens it- 
self will obscure any background objects directly behind it. 
In addition, the flexion profile described by equation [l9] be- 
comes very large as a; ^ 0. A more reasonable model of the 
lens profile might be a softened isothermal sphere, in which 
the convergence is given by 



k{x) = 1 + 



29% 



1 + 



(21) 
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Figure 3. Top Panel: The expected signal to noise for the filter 
described in equation [m] plotted as a function of fi and R. Bot- 
tom Panel: The expected signal to noise for the family of filters 
described in equation |16| as a function of I and R. 

which gives rise to a flexion signal given by 

Figure |3] shows the expected signal to noise for this soft- 
ened isothermal profile for the filter functions described by 
equations [14] and |16[ For each, we assume a velocity dis- 
persion of (T„ — 500km s~^, which gives rise to 9e — 5". 95 
assuming Dig/Ds = 0.8. We take n = 35 arcmin"^ and 
ajr = 0.1/". This is a low value compared to the disper- 
sion in flexion measurements found by Goldberg and Bacon 
(2005) or Goldberg and Leonard (2007); however the simu- 
lated data with which we are concerned in this paper has a 
lower dispersion than that seen in real data, thus the value 
of ajr taken here is representative of our data. In addition, 
the measured values found in Goldberg and Bacon (2005) 
and Goldberg and Leonard (2007) are likely to be overes- 
timates of the "true" flexion dispersion, as these measure- 
ments were carried out in lens flelds, thus we expect to see 



a flexion signal in these background galaxies. We flnd the 
peak attainable signal to noise for each fllter to be ~ 18, 
with a slightly higher peak signal to noise found using the 
piecewise-continuous flexion filter function. 

It is important to note that the expected signal to noise 
ratio has been computed here under the assumption that 
there is only one lens present in the regime being considered. 
However, clusters of galaxies tend to be quite clumpy; thus 
we expect the signal to noise achievable in clusters to be 
somewhat lower, though still appreciable. 

It is clear from the figure that the signal to noise prop- 
erties of these two fllter functions differ signiflcantly. Thus 
comparisons between the aperture mass measured using dif- 
ferent fllters provides a robust method for checking the ac- 
curacy of these maps. It is clear that for most values of ui, 
the piecewise-continuous function favours smaller apertures; 
that is, apertures more closely matched with the character- 
istic size of the lens. The polynomial fllters, on the other 
hand, favour larger apertures as the polynomial order is in- 
creased, and indeed should provide robust detections for a 
broad range of aperture sizes for any given choice of poly- 
nomial order I, particularly if I is large. 

In addition, for a given choice of aperture size, increas- 
ing I should decrease the minimum scale on which struc- 
tures are resolved, as increasing / decreases the radius at 
which the fllter function reaches its peak. Similarly, decreas- 
ing the value of vi for the piecewise-continuous filter should 
decrease the minimum scale on which structures can be re- 
solved. Thus it is possible with both these filters to resolve 
structures on both large and small scales by appropriately 
tuning the filter. This is important in clusters of galaxies, 
where one is interested in both the large-scale structure of 
the cluster potential, but also the smaller-scale substructures 
within it. 



4 APPLICATION TO SIMULATIONS 

We test our method on two difi'erent simulated lens systems 
using a ray-tracing simulation to artificially lens background 
test galaxies through a given convergence field. 



4.1 Ray-tracing Method 

To carry out the ray-tracing, we first consider the mapping 
between the source and lens planes. This is given by /3 = 
— a{0), where /3 are the coordinates in the source plane, 
are the coordinates in the lens plane, and a are the deflection 
angles evaluated as a function of the lens plane coordinates. 
In Fourier space, the convergence and the deflection angle 
are related by 

= - t.2^f p ^- (23) 

Thus, in order to calculate the deflection angles arising 
due to a given lens convergence, one simply needs to Fourier 
transform the convergence fleld, apply equation |23| and then 
compute the inverse Fourier transform. In order to achieve 
sufficient numerical accuracy and avoid edge effects, zero- 
padding is applied surrounding the convergence field. For all 
the simulations described in this section, we use a padding 
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or space-based observations. Further, current and upcoming 
surveys and instruments will regularly be able to achieve a 
background source count of in excess of 35/arcmin^. Some 
examples are the Large Synoptic Survey Telescope (LSSlQ, 
which should offer 50 sources/arcmin'^, the Canada- France- 
Hawaii Telescope (CFHT) Legacy Survey, which obtains 35 
sources/arcmin^ (Gavazzi & Soucail 2007), and Suprime 
Cam, for which conservative estimates quote 35/arcmin^ 
and higher (van Waerbeke et al. 2006). 

The pixel locations for each source are drawn from a 
uniform random distribution, while each of the two compo- 
nents of the complex ellipticity, e^, are drawn from a random 
Gaussian distribution with a mean of zero and a standard 
deviation of 0.2. The background sources are all placed at 
the same redshift (zs — 1.0), and the brightness in the source 
plane generated according to 

/(a;,y) = ^e-K^^)\ (24) 

i 

where the sum is taken over all background galaxies, 

J ((1 - e^)x' - e2y'f + i-e2X' + (1 + ei)y')' ^ , 
n = '(25) 

and x' = X ~ Xi. We take rgai = 6kpc to be a characteristic 
scale for the galaxies, and define a — rgai/1.22dgai, where 
Da is the angular diameter distance to the source plane. 

The background galaxies are lensed by computing the 
mapping from each pixel location in the lens plane to the 
source plane and using a cubic interpolation to determine 
the brightness distribution at that location. 



4.5 Computing the Aperture Mass Statistic 

To compute the aperture mass over the field of view, we 
use a grid of 1000 x 1000 apertures of a fixed radius. The 
flexion is measured using the weighted HOLICs formalism 
laid out in Goldberg & Leonard (2007). As in Leonard et 
al. (2007), we reject any flexion measurements for which 
a\^\ > 0.2, where a is the measured semi-major axis of the 
lensed galaxy. This is done to ensure that no strongly lensed 
objects are included, and to reduce contamination from bad 
measurements or blended sources. 

In order to compute the noise in the aperture mass mea- 
surement, the data are randomised by rotating each flexion 
vector by an angle drawn randomly from a uniform distri- 
bution in the range ^ 9 ^ 2n and computing the aperture 
mass on the randomised data. This procedure is carried out 
5000 times, and the noise is given by the standard devia- 
tion of the aperture mass maps generated using randomised 
data. 



5 RESULTS 

5.1 Single Lens Galaxy 

As a first test of the aperture mass method, we consider here 
a single lens galaxy with a velocity dispersion of 500km s~^, 

1 LSST Dark Energy Task Force Whitepaper; 



on all four sides of the image that is 3.5 x the size of the 
input image. 

4.2 Analytic Simulation 

As a first test of the aperture mass technique, a purely ana- 
lytic lens convergence is generated by laying down a number 
of randomly distributed, circular test "galaxies" with a con- 
vergence defined by equation |21| and with a pixel scale of 
0".5 pix~^. For these simulations, we take Dis/Ds — 0.8, 
thus for a source redshift of Zs — 1.0, the lens will be at a 
redshift of zi ~ 0.16. 

We test our method on both a single lens galaxy and a 
group of galaxies. The velocity dispersion for the single lens 
case is taken to be a„ = 500km s~^, while for the group 
simulation, the velocity dispersion of each lens galaxy is 
drawn randomly from a Gaussian distribution with a mean 
of 220 km s^^ and a standard deviation of 20kms^^. The 
pixel locations of each test galaxy in the group are also 
drawn from a Gaussian distribution centred on the centre of 
the image grid, and with a standard deviation set to A^pix/10, 
where A'pix is the width of the image in pixels. 

4.3 N-body Simulation 

The advantage of using a purely analytic potential is that the 
method can be tested and checked for accuracy at various 
stages of the analysis pipeline, as the expected flexion and 
deflection angles can be calculated analytically. However, 
this method does not generate a particularly realistic model 
cluster of galaxies. 

For this purpose, we make use of a cluster extracted 
from the Millennium simulation (Springel et al. 2005) at a 
redshift of zi = 0.21. The resolution of the cluster image 
used is 5h^^ kpc/pixel (2". 07 pix~^), corresponding to the 
gravitational softening of the simulation, and the projections 
are computed using a nearest grid point interpolation. The 
cluster has a characteristic radius -R200 = 2.0 h~^Mpc and 
a mass of A/200 = 1-23 x 10^^ H'^Mq. We extract the cen- 
tral 1.0/i^^Mpc of this cluster for use with our ray-tracing 
simulation. 

It is important to note that, when carrying out the ray- 
tracing using this cluster as our lens, we do not apply any 
smoothing or inter-pixel interpolation prior to computing 
the deflection angles across the lens plane grid. As the cluster 
data is quite noisy, and the pixel scale is significantly larger 
than that used for the simulated background galaxies, we 
expect this inter-pixel noise to result in noisier fiexion data 
than that found when using a lens simulated as described in 




4.4 Background Galaxy Population 

For all the ray-tracing simulations, we use an image resolu- 
tion of 0.l"pix~^, and a background number density of 35 
sources/arcmin^. This may seem to be a rather optimistic 
value compared to typical ground based observations, and 
indeed when compared to the analysis of Abell 1689 carried 
out by Okura et al. (2008) using Subaru data, in which they 
used 5 sources/arcmin^. However, the Leonard et al. (2007) 
analysis used 75/arcmin'^, which is rather high, but typical 
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Figure 4. The flexion aperture mass signal to noise contours for a 
simple softened isothermal lens and randomly-sampled theoretical 
data overlaid on a density plot of the underlying convergence. 
The solid contours show positive signal to noise, while the dashed 
contours show regions of negative signal. The zero point contour 
is clearly labelled. 

corresponding to an Einstein radius oi 8e = 5". 95 assuming 
Dis/Ds = 0.8. In addition, we consider data in which the 
flexion is measured exactly; i.e. we randomly sample points 
on the lens plane grid at a rate of 35 arcmin"'^, analytically 
compute the flexion at these points, and use these data to 
compute the aperture mass. Any data points that are found 
to lie within the Einstein radius of the lens are discarded 
prior to computing the aperture mass. 5000 randomisations 
of the data are carried out, the resulting dispersion giving a 
measurement of the noise associated with the aperture mass. 
We use the polynomial flexion fllter function described in 
§ |3.1.2[ taking 1 = 3 and R = 60". 

Figure |4] shows the signal to noise map for the aperture 
mass computed in this way. As expected from the calcula- 
tion in § |3.2[ we flnd the peak signal to noise for this recon- 
struction to be Speak = 18.1. A very noticeable feature of 
the aperture mass map generated in this way is that there 
is a clear radius at which the signal becomes negative. We 
expect the aperture mass to become negative at large radii 
as the mass filter function becomes negative, and the radius 
at which this occurs will be dependent on the mass profile 
of the lens, the strength of the lens, the shape of the fil- 
ter function used, and the size of the aperture. This implies 
that if this zero contour can be resolved in real images, and 
a reasonable shape for the mass profile of the lens can be 
assumed, one might be able to infer the mass of the lens (or 
its Einstein radius) simply by measuring the radius of the 
aperture mass zero contour. 

A detailed treatment of this idea is beyond the scope of 
this paper; however, as a test of this concept, we measured 
the zero point radius of the aperture mass for a simple sin- 
gle softened isothermal lens of variable Einstein radius, and 
for various combinations of aperture scale and polynomial 
order. Figure |5] shows the zero point radius plotted as a 
function of 9e/R for an aperture radius of 7? = 120" and 
a polynomial filter with i = 5. It can clearly be seen from 
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Figure 5. The zero point radius iJo as a function of O^/R for 
an aperture size of R = 120" and a polynomial order I = 5. The 
discrete points represent the actual data, and the solid line shows 
the line of best fit for the data. The solid curves represent regions 
of positive signal, while the dashed curves show areas of negative 
signal. 

the figure that the relationship between i?o and 9e appears 
to be linear. Different combinations of / and R yield similar 
linear results, though with slightly different slopes as I is 
varied. 

Clearly as the lens model becomes more complicated, 
or less smooth, the relationship between the Einstein radius 
and the zero point radius will become less straightforward. 
However, it does seem that if these zero point contours can 
be resolved using the flexion aperture mass, this might pro- 
vide a straightforward way to put at the very least broad 
constraints on the mass contained within the structures seen. 

5.2 Analytic Simulation of a Group of Galaxies 

We now consider a group of softened isothermal profile 
galaxies generated as described in § |4.2| These galaxies were 
generated to lie within a region corresponding to a physical 
scale of 0.67Mpc on a side at a redshift of zi = 0.16. The 
aperture mass statistic was calculated using both the poly- 
nomial and piecewise continuous filter functions described 
in § |3.1.2| and § |3.1.1| In addition, our flexion data are gen- 
erated using the method described in § |4.4[ 

Figure [6] shows the aperture mass signal to noise con- 
tours computed using the piecewise-continuous filter func- 
tions of § |3.1.l"] using an aperture scale of i? = 90" and taking 
I'l = 0.05 and 0.1. With both values of the inner radius, the 
signal to noise contours are seen to quite accurately trace 
the shape of the convergence, to a peak signal to noise of 
Speak = 5.0 for vi = 0.05, and Speak = 5.1 for u\_ = 0.1. 
As expected, the signal to noise contours for the smaller 
inner radius resolve structures on smaller scales than that 
for larger inner radius. Indeed, reconstructions involving a 
different aperture size, though not shown here, exhibit the 
expected behaviour in terms of the scale on which structures 
are resolved. In addition, a zero point contour is quite clearly 
resolved in both signal to noise plots. 

Figure [7] shows results obtained using the polyno- 
mial filter functions, shown for the combinations [R, I] = 
[60", 3], [90", 5] and [90", 10]. The peak signal to noise ob- 
tained in these reconstructions is Speak = 4.1, 4.5 and 4.2, 
respectively. Here again the signal to noise contours accu- 
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Figure 6. The aperture mass signal to noise contours computed 
using simulated data and overlaid on an image of the under- 
lying convergence. The aperture mass is computed using piece- 
wise continuous filter functions with ui = 0.05 (Top Panel) and 
ui = 0.1 (Bottom Panel). For each, the outer radius is taken to 
be U2 = 0.95. 



rately trace the shape of the galaxy group, and are able to 
resolve structures on subarcminute scales. In addition, as 
expected, we see that increasing the aperture size increases 
the minimum scale on which substructure can be resolved, 
and increasing the polynomial order allows us to probe these 
smaller-scale substructures. Further, we again resolve a zero- 
point contour, the radius of which, when compared to the 
linear relation shown in § |5.1[ yields an Einstein radius which 
is broadly consistent (within a factor of ~ 2 — 3) with that 
found by modelling all the mass contained in the lens as a 
single softened isothermal sphere. 

As a check on the performance of this method, we com- 
pare our results to a weak lensing mass reconstruction, by 
using the measured ellipticities of our lensed images to per- 
form a finite field inversion following the method of Seitz 
and Schneider (2001). The reconstruction was performed on 
a 100 X 100 grid, with a smoothing scale of 120 pixels (0.048 
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Figure 7. The aperture mass signal to noise contours computed 
using simulated data and overlaid on an image of the underlying 
convergence. The aperture mass is computed using polynomial 
filter functions with 1 = 3, R = 60" (Top Panel), 1 = 5, R = 90" 
(Middle Panel), and I = W R = 90" (Bottom Panel). The solid 
curves represent regions of positive signal, while the dashed curves 
represent regions of negative signal. 
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Figure 8. A reconstruction of the convergence from shear data 
using the finite field inversion technique of Seitz & Schneider 
(2001); see text for details. 

times the image size) - chosen to be large enough to include 
averaging over sufficient background galaxies to lessen the 
noise from their intrinsic ellipticity dispersion. The resulting 
mass reconstruction is shown in figure|8] This reconstruction 
is clearly unable to resolve the small-scale structure of this 
galaxy group, or to accurately trace the overall shape of the 
convergence distribution. 

5.3 Cluster of Galaxies from N-body Simulations 

Finally, we consider a cluster of galaxies taken from the Mil- 
lennium simulation and described in § |4.3| Figure [9] shows 
the signal to noise contours for the flexion aperture mass 
statistic applied to this cluster using the polynomial filter 
functions with « = 5 and i? = 60", 90" and 120". The peak 
signal to noise for these maps is found to be Speak = 3.8, 3.0, 
and 3.3, respectively. Similar results were obtained using the 
piecewise-continuous filter functions. 

We expect, and indeed find, the signal to noise map 
to be somewhat noisier than that seen with the analytic 
simulations. This manifests as spurious detections on small 
scales of structures that do not exist, which undoubtedly 
result from using rather noisy data to compute the deflection 
angles. As expected, when the aperture scale is increased, 
these detections are removed, as the aperture mass becomes 
less sensitive to structures on those scales. 

While considerably noisier, the flexion aperture mass 
does still accurately trace the shape of the central part of 
the cluster, and resolves two smaller substructures (one to 
the left and one directly below the centre of the cluster) in 
all three signal to noise maps. However, as found in Leonard 
et al. (2007), the flexion measurements are generally quite 
insensitive to the smooth component of the cluster potential, 
so outside of the central, cuspy region, the flexion signal 
drops off rather quickly. 

Again, we can compare the results from flexion with 
that from shear. As for the analytic group, we use the Seitz 
and Schneider (2001) finite field inversion technique to re- 
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Figure 9. Aperture mass signal to noise contours computed us- 
ing simulated data lensed through a cluster from the Millennium 
simulation, overlaid on the underlying convergence. Only positive 
signal to noise contours are plotted. The aperture mass was cal- 
culated using a polynomial filter function with I = 5, and with 
aperture radii of i? = 60" {top panel), R = 90" {middle panel), 
and R = 120" {bottom panel). 
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Figure 10. Shear reconstruction of the convergence field using 
the finite field inversion technique described in Seitz and Schnei- 
der (2001); see text for details. 

construct the convergence from our shear data. The mass 
was reconstructed on a 100 x 100 grid, with a smoothing 
scale of 0.029 times the image size. The mass reconstruction 
is shown in figure |10| In this case, the shear does pick up 
the overall shape of the cluster quite accurately, but does not 
clearly resolve the smaller substructures seen in the flexion 
aperture mass signal to noise map. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we have derived an aperture mass statistic for 
(first) flexion, and shown that it provides a robust method 
for detecting structures within clusters of galaxies on sub- 
arcminute scales. Whilst the filter functions used here are by 
no means optimal, they offer an excellent proof of concept 
for this technique, and show that even using non-optimal fil- 
ters, and a moderate number density of background galaxies, 
mass concentrations can be detected on a range of different 
scales with an appreciable peak signal to noise that varies 
very little with changes to the aperture scale radius or filter 
parameters {vi, V2, and I). 

We have also shown that the zero point contours of the 
aperture mass map can be clearly resolved in simulated data, 
and described briefly how these zero points might be used 
to provide an estimate of the mass contained in structures 
identified using this method. This is, of course, model- and 
filter-dependent, but might provide a good starting point for 
maximum-likelihood modelling in clusters of galaxies where 
fiexion constraints are incorporated to describe the substruc- 
ture within the cluster, for example. 

It is clear that the filters used in the preceding analysis 
are not optimised for any particular lens, but can be used 
to detect structures with a range of different mass profiles. 
This will be useful for detecting substructures where one has 
little or no prior knowledge of the mass profiles. However, 
it is clear that improvements can be made to the method, 
and improved signal to noise obtained, by using filters that 
are tailored to the structures being studied. This will be 



particularly important when dealing with noisier data, or 
data for which the background source number density is low, 
thus reducing the maximum attainable signal to noise. 

Our understanding of the physical properties of galaxy 
clusters and their utility as cosmological probes will be 
greatly enhanced with the advent of large surveys such as the 
Dark Energy Survey (DES), augmented with complemen- 
tary data from experiments such as the South Pole Tele- 
scope (SPT). For the large cluster samples, gravitational 
fiexion will provide constraints on the mass substructure of 
clusters at a higher resolution than weak lensing, and probes 
greater distances from the cluster centre than strong lens- 
ing. Clearly, a priority is to simultaneously use lensing data 
from various regimes in conjunction with other types of data 
to obtain a detailed multi-scale view of the largest bound 
structures in the universe. 
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APPENDIX A: RELATIONSHIP BETWEEN 
THE CONVERGENCE AND FLEXION 

The convergence can be related to the measured flexion in 
Fourier space by (Bacon et al. 2006): 
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iko 
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T2 



(Al) 



This implies that the real-space relationship between k 
and the flexion is a convolution: 



.(X) = 



y d^x'j^(x')s>(x - x') 



where Er is a convolution kernel to be determined. 
Noting that 
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where _F(/(k)) denotes the fourier transform of the function 
/(k), we find that 

Et{X) = ^, (A5) 

where X — x\ ^ 1x2 = xe*'*, and is the phase of the 
complex number X. 



APPENDIX B: EXPRESSING THE APERTURE 
MASS STATISTIC IN TERMS OF THE E-MODE 
FLEXION 

The aperture mass is defined, in terms of the measured flex- 
ion, as 



m(xo) : 
5R 



y d^x'jr(x') y d^xS>(x-(-Xo-x')w(|xj) 



(Bl) 



Making the transformation y = x' — xq, and writing d'^x — 
X dx d<j), we obtain 



m(xo) ^ 
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where X — xi + ix2 — \x.\e"^. We consider the rightmost 
integral first, which can be rewritten as 

r2vr J J, 
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This integral has a regular singularity a,t X = Y , and us- 
ing the transformation dX = iXd()), can be expressed as a 
contour integral: 

i r dx 
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2n J X(X-y)' 
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where the integration is carried out along a circular contour 
of radius This integrand has two simple poles at X = 
and X — Y. The residues at these singularities are given by 
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Thus, applying the residue theorem, 
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where 
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Returning to equation |B2| we find that 
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Now, the E-mode flexion is given by 



thus we can write 



m(xo) = - / d'^y Tsiy - :x-o) - 
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which leads directly to equations [9] and |10[ 

Thus, we have expressed Qr(x) in terms of 10^^(2;) as 



y w{y) dy. 



(B12) 



We now aim to derive the reverse relation; i.e. to express 
wjr{x) as a function of Qjr{x). 

Recall that first flexion is related to the convergence via 



T — Ti + iT2 — Ok — {di + id2)K. 
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It is useful to transform the derivative operators into polar 
coordinates as 



di — cos 



02 = sm0c*r H 04, 



(B14) 



Using the above expressions, and recalling that the E-mode 
flexion is defined as J^e = ^\ cos (f>+J-2 sin (j), we flnd that the 
E-mode flexion is related quite simply to the convergence via 
J-E = drii- Thus, we can express the aperture mass statistic 
as 



i(xo) = /d'x(|^ j Q, 
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r— dr. 
or 
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Integrating the inner expression by parts and simplifying, 
we find that 

m(xo) = ~ Jd'^ ^(^ + xo) (^lQr{r) + . (B16) 

Comparing this expression to equation |4] it follows that 

^Ar)-~lQAr)-'^- (B17) 



